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In this paper, we study the experimental signatures of a gluophilic Z' at the LHC, in 
particular through the analysis of three jets or four tops events. The Z' couples to gluons 
through dimension six operators and the parameter space is constrained with experimental 
searches released at 7 and 8 TeV by CMS along these two different channels. Existing 
constraints coming from the study of dark matter where the Z' represents a possible mediator 
between the latter and the Standard Model are also included for comparison. Prospects at 
= 13 TeV allow us to evaluate for which values of the parameter space a gluophilic Z' 
could be discovered during the next run of the LHC. In particular, we show that the analysis 
of the three jets invariant mass could provide a clear signal (> Str) for masses of the Z' 
above 300 GeV. Four tops events bring in addition further discovery potential for heavy Z' 
(above ^ 2 TeV). A combination of both signals in four top channels and three jets analyses 
during the next run of the LHC could thus provide a clear signal of the presence of a heavy 
gluophilic Z'. 
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I. INTRODUCTION 

One of the most simple extension of the Standard Model (SM) is to provide the latter an 
additional U{1) abelian symmetry [T]. The associated gauge boson - usually denoted by Z' in the 
literature - has been given a particular attention in the last decades in particular as a potential 
candidate for mediating interaction between the dark sector and our visible world |2|. Within this 
approach, the Standard Model fermions can be considered to be charged or not under the additional 
- so called U{l)x in this paper ~ gauge group. In the case of charged fermions, a particular care 
must be devoted to anomaly cancellation and flavour changing constraints. The B — L models are 
among the most popular example, and satisfy these requirements by considering a very heavy Z'. 
On the other hand, string inspired models propose an anomaly cancellation a la Green-Schwarz (SHS] 
where a very light U{l)x mediator is generated with the use of axionic couplings and generalized 
Chern-Simons terms, or via Stueckelberg realization of the Z' boson |7|. Another option is to 
consider neutral Standard Model fermions under the new gauge group U{l)x [TO] . 

By adding heavy degrees of freedom, charged under both SM and U{l)x symmetries, one 
can generate effective operators of dimension six and compute an effective interaction Lagrangian 
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between a Z' boson and the SM gauge bosons. Such scenarios have already been studied in 
the context of Dark Matter (DM) model building |ini413j . In such framework, dimension six 
operators are suppressed by a factor (M being the mass of heavy fermions integrated out) 
and the dark matter sector, charged under U{l)x, couples with the weak or coloured sector of the 
Standard Model. An interaction with the weak SM sector was shown to have possible astrophysical 
signatures mm, while a coupling to the coloured SM sector was constrained m using mono¬ 
jets events at the LHC m, as well as indirect detection constraints arising from astrophysical 
measurements. 


In this paper we focus on the case aforementioned where Standard Model fermions are neutral 
under U{l)x and where the Z' boson interacts with the Standard Model gluons through effective 
operators of dimension six (and possibly to a dark matter sector). Interactions between the Z' 
mediator and the SM gluons will be detailed in Sectionjllj The presence of such a feeble coupling has 
interesting features in the quark physics which may be discoverable during the Run-2 (^/s = 13 TeV) 
of the LHC. Therefore, this study aims to show what sensitivity can be reached for the theoretical 
model developed in m given the latest LHC Run-1 [y/s = 7 or 8 TeV) experimental results. 
For completeness, the dark matter constraints and prospect studies for the LHC Run-2 are also 
included. In Section we discuss the possible experimental signatures. We notably identify two 
interesting and complementary channels that are analyzed in more details in Sections IV and |Vj 
where we use the existing experimental constraints to investigate what would be the maximal 


coupling allowed for the Z' interaction with the SM gluons. Finally, in Section VI we present our 
conclusions concerning the potential of discovery of the model in the next years. 


II. THE GLUOPHILIC Z' MODEL 

As mentioned in the introduction, a new U{l)x gauge group is added to the SM under which 
SM fermions are considered to be neutral [TO]. Effective operators of dimension six between the 
Z' and the SM gluons are then assumed to be generated at the loop level by integration of heavy 
fermions - namely ^l,r of mass M - charged both under U{l)x and SU{3)c gauge symmetries. 
Such loops produce at low energies effective interactions between the Z' bosons and the SM gluons 
as depicted in Fig.[^ The heavy fermions 'kget mass through a spontaneous symmetry breaking 
mechanism of the U{l)x symmetry with the use of a heavy Higgs boson field (<k). A gauge invariant 
Lagrangian describing such a theory can be written by realizing the gauge symmetry non linearly, 
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a la Stueckleberg, as shown in what follows. 




G 
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FIG. 1: .^’coupling to colored sector of the Standard Model. The heavy mediators are suppressed when integrating 
over all considered heavy states | 10 |. 


The heavy Higgs field can be written a la Stueckelberg 
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if the Higgs boson mass is assumed to be much heavier than the Z' boson. Only the axionic 
component of the initial field thus remains in the low energy theory. The Z' and the axion field 
transform non-linearly under [/(l}x as follows 


SZ'^ = dfj,a 


SGx = ^ a where 9x = ■ 

^ V 


( 2 ) 


Moreover, in order to write down a gauge invariant interaction lagrangian, the covariant derivatives 
of the axion Ox and SM gluon fields are introduced 


DnOx = dnOx - ^Z' 


^ 2 


abc/^b /oc 


= d.G^g + gr'^'^GlG: 


OL^ 




(3) 


where gx denotes the U{l)x gauge coupling and g stands for strong coupling constant. The U{l)x 
invariant Lagrangian under Equation Q can finally be written as follows 


^ = ^SM + ^ dgd^DgOxTriGG) + d'gd^D-exi:T{G^pGP)+ 


e'DPOxTviG.pVgGPn + egT»^0xTr(G«.P"G^“) 


(4) 


where Csm is the Lagrangian describing the Standard Model interactions, and the dual field- 
strength GP'^ = dZ^Gpa is introduced to protect the CP parity of the previous operators. An 
explicit computation of fermionic loops has been released in |10] where it has been shown that. 


interestingly, only the operators proportional to dg and Cg are actually present in the theory and 
related as follows 


Cg — 2dg. 


( 5 ) 
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Z' 300 GeV 

Z' 500 GeV 

Z' 800 GeV 

Z' 1.5 TeV 

Z' 3 TeV 

Z' ^ qqg 

2.0 MeV 

25 MeV 

0.25 GeV 

5.3 GeV 

0.16 TeV 

Z' —y ttg 

- 

0.27 MeV 

19 MeV 

0.82 GeV 

30 GeV 

Z' —)■ ttti 

- 

- 

33 eV 

0.79 MeV 

0.13 GeV 


TABLE I: Decay widths of several Z' decay channels, for various Z' masses and an effective coupling 
dg/M'^ = 10“®. The branching ratio in the multijets channel Z' —> qqg always dominates almost exclusively. 


Thus, only these two interaction terms will be considered in this stndy, the only free parameters 
remaining in the model being dg/M‘^ and the coupling constant gx- Finally, for the Z'gg interaction 
the interaction terms presented in Equation Q can be written explicitly 


C D 

D 


dg 

w 

dg 


M2 


df^D^OxTriGG) - 2M^D^0xTi{Go,,V’'&^) 
gxd^Z'^e^^-^P-d^G^dpG^ - gxZy^P^d^^Gi^d^dpG^ 


( 6 ) 


and the two associated vertex functions for processes involving Z'{pz')G{pi)G{p 2 ), symmetrized 
with respect to the two gluon fnnctions, are simply 








+ 




+ (P2)m(pi)"(pi)re^™’'" 


(7) 

(Pl-P2)(p2-Pl)re^""" 

( 8 ) 


where Uout is the number of outgoing particles in the process considered. The gauge coupling gx 
can generically take values of order 0(0.1 — 1). For simplicity, gx will be fixed in what follows to 
unity in order to incorporate it into the definition of the free coupling of the model, dg/M"^. 

As far as constraints coming from dark matter are concerned, we will see in the next sections 
that the study of [TO] provides stringent constraints on the conpling dg/M'^ for low masses of the 
Z' nsing an analysis of the monojets + missing Et at the LHC (8 TeV). Such constraints will be 
only mentioned in the following study as a possible restriction of the resnlts, since the presence 
of a gluophilic Z' does not necessarily implies the presence of a dark sector - nor constitutes the 
exclusive possibility of interaction of DM with the SM. 


III. PHENOMENOLOGY AT HADRON COLLIDERS 


Thanks to its coupling to the colored sector, the Z' boson phenomenology at hadron colliders for 
these gluophilic models is pretty rich and in some aspects distinct from other BSM physics scenarios. 
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FIG. 2: Example of Feynman diagrams for the processes involving a Z' boson discussed here. 

Z' bosons may be produced directly in proton-proton collisions, or enhance the Standard Model 
cross-sections of rare processes through offshell contributions. The leading order process gg Z' 
is suppressed due to the Landau-Yang theorem and the typically narrow Z' width (cf Table ; 
therefore direct production may occur through the loop-induced process qq^ Z' ^ or in association 
with a parton from initial state radiation (ISR), qg —)• qZ' and gg —)• gZ'. In this study we focused 
on the associate Z' + ISR jet production, the related cross-section being easier to evaluate. 

The decay of a tree level Z' boson into two gluons is forbidden for the same reason as the 
gg —)• Z’ process is. Possible decay channels thus involve final states with at least three quarks 
or gluons (e.g. Z' —)• qqg), including the experimentally interesting particular case of top quarks, 
and possibly radiated electroweak or Higgs bosons. For a more quantitative understanding we 
determined the Z' decay widths in various channels, a few of which are summarized in Table 
For that we relied on the FeynRules 2.3 package m to establish the Feynman rules corresponding 
to the effective Lagrangian Q, allowing the computation of leading order matrix elements and 
cross-sections by the MadGraph 5.2.2.3 Monte Carlo generator m- One can observe that the 
channel Z' —)• qqg is largely dominant, the branching ratios of channels involving top quarks or W 
bosons being typically below 1% apart for large Z' masses. While the latter provide resonant final 
states involving top quarks or electroweak bosons which are very clean experimental signatures, we 


^ Note that no loop processes have been included here. Ideally, a full QCD computation including virtual corrections 
would be required for a better precision on the results. 
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Z' 300 GeV 

Z' 500 GeV 

Z' 800 GeV 

Z' 1.6 TeV 

Z' 3 TeV 

= 7 TeV 

58 pb 

16 pb 

3.6 pb 

0.14 pb 

0.56 fb 

Vs = 13 TeV 

0.67 gh 

0.22 

69 pb 

7.8 pb 

0.38 pb 


TABLE II: Cross-sections otpp —>■ Z'j associate production for different Z' masses, with an effective coupling 
dg/M^ = 10 - 6 . 

focused on the (light flavored) multijets signature, which is very competitive thanks to its branching 
ratio close to unity. 

We also studied non-resonant Z'-mediated contributions to rare SM processes, which may pro¬ 
vide a nice complementarity to the direct production, in particular to probe very large Z' masses. 
One interesting such example is the four-tops production pp —)• titi, which can be significantly 
enhanced by these new contributions. 


IV. SENSITIVITY IN THE MULTIJETS CHANNEL 

We focus in this section on the experimental signature corresponding to the associate production 
of a Z' boson and an ISR parton, with a Z' assumed to decay into qqg with a 100% branching 
ratio. Various searches by the ATLAS and CMS collaborations look for the resonant production 
of heavy particles decaying into multijets final states [MSI]. However most of them are irrelevant 
here since they analyze either dijet production, or assume pair-production of the heavy particles 
(e.g. pp —)• gg with a gluino RPV decay). More useful are the measurements of the differential 
cross-section of QCD three-jets production at ^/s = 7 TeV, as a function of the three-jet invariant 
mass |23l I24| . We evaluated the sensitivity of the measurement by CMS [23] to a potential Z' 
signal, and reinterpreted the measured cross-section into upper limits on the Z' effective coupling 
to gluons dg/M‘^. 

The CMS measurement selects events with at least three jets {pp > 100 GeV, \y\ < 3.0) and 
provides the observed cross-section as function of the invariant mass rrijjj of the three leading pp 
jets in the range 445 < rnjjj < 3270 GeV, in two bins of rapidity |y|max < 1 and 1 < |y|max < 2, 
where |?/|max corresponds to the largest rapidity among the three considered jets. To evaluate the 
Z' signal acceptance for this selection and construct the corresponding three-jet mass spectrum, 
we generated Monte-Carlo samples of signal events for different Z' masses (and a decay width fixed 








8 




(b) Differential fiducial acceptance for pp —> Z'j —>■ qqgj 


FIG. 3: Three-jet invariant mass spectrum for QCD (left) and Z' signal models of various masses (right). 


to^ 1 GeV) at leading order using MadGraph 5.2.2.3[I7] with the CTEQ6L1 set of parton distri¬ 
bution functions [25], and interfaced to Pythia 6.4.28 [26] for parton showering, hadronization and 
modeling of the underlying event. Jets were then reconstructed from visible particles (no detector 
simulation involved), with their momentum randomly smeared by 10% to mimic the finite detector 
resolution. Apart from this emulation of the jet energy resolution, the detector reconstruction and 
data acquisition inefficiencies were neglected - a reasonable assumption as they are at the percent 
level [23] [27]. The signal production cross-sections were evaluated altogether, a few typical val¬ 
ues are gathered in Table [^ The fiducial acceptance corresponding to the aforementioned event 
selection was found to increase from 5 to 50% for Z' masses in the 300 — 3000 GeV range. The low 
acceptance at small masses originates from the rather hard jet pp cut used in the selection. 

A quick look at the invariant mass spectrum for the signal, shown on Fig. |3(b)] already provides 
some useful insight. While one would a priori expect that the presence of a resonant peak^ at the 
Z' mass on top of an exponentially decreasing QGD background would provide the best evidence, 
the presence of the ISR parton complicates the situation. Indeed, it tends to be very energetic 
and be selected as one of the jets entering the invariant mass computation. As a result the signal 
invariant mass spectrum is not a sharp peak, but considerably smeared, even more so for light Z'. 
This is in fact a rather nice feature, as it allows some signal to populate bins at large invariant 
mass (above the resonance mass) where the QCD background is much lower. On the other hand, 

^ Note that for such values of the decay width - as computed in Table |^- the influence of the latter on the simulation 
turns out to be neglectable. 

® The Z' narrow width and the good jet energy resolution in ATLAS or CMS would allow such a feature. 
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FIG. 4: Current LHC sensitivity to the .Z'-gluon coupling in the three-jets and four-tops channels, together 
with dark matter-related constaints from m- 

it would make the interpretation of a hypothetical observed excess less straightforward. One 
can see that the jets tend to be more central in signal events than in QCD events, therefore the 
measurement in the bin |y|max < 1 can be expected to provide most of the sensitivity. 


Exclusion limits on Z' signal 

As the CMS measurement showed good consistency between the observed data {^/s = 7 TeV, 
L = 5.0 fb“^) and the Standard Model QCD predictions at NLO -|- non-perturbative corrections 
(Fig. |3(a)D , we reinterpreted these observations in terms of upper limits on the Z' free coupling 
to gluons dg/M"^, as a function of the Z' mass, using the signal three-jet invariant mass spectra 
determined as described in the previous paragraph. For that we relied on the public 2.0 version of 
the HistFitter software ^ [28], which allows the computation of exclusion limits at 95% conhdence 
level in the CLg formalism with a test statistic built from a one-sided profile likelihood ratio. The 
limits were determined by a simultaneous maximum-likelihood estimation (MLE) fit of the observed 
invariant mass distributions in the two |y|max bins with the signal and background components, 
with a free signal strength. We observed that the contribution from the bin 1 < |y|max < 2 is 
indeed minor, improving the upper limit on the signal cross-section by only 10%. The uncertainty 


This software is largely used by the ATLAS collaboration, notably for all SUSY searches results. 
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on the background prediction included in the likelihood is detailed later on, and amounts to 20% in 
the most sensitive bins. Finally, the tool was setup to use asymptotic formulas for the distribution 
of the test statistic, as explained in [25]. 

The upper limits thus determined are shown on Fig. together with the existing dark 
matter-related constraints established in m, which set limits on the parameter dg/M"^ based 
on searches for events with monojets + missing Et at LHC as well as astrophysical constraints 
coming from indirect detection and relic density. As a matter of fact, both approaches reach 
sensitivity to free Z'-gluons couplings of the same order of magnitude, the reinterpretation of the 
CMS measurement proposed here performing a bit better at high Z' mass (above 1 TeV) while 
the dark matter search dominates at lower masses. 

Projections for the next LHC run at ^/s = 13TeV 

The next phase of data-taking at the LHC (2015-2018) will provide a new dataset with increased 
luminosity (up to 100 fb“^) and a center-of-mass energy raised to 13 TeV, which will particularly 
improve the sensitivity to heavy Z'. We evaluated the discovery potential that could be achieved 
during this exploitation phase. To predict the QCD background yields in these new conditions, 
we generated Monte-Carlo samples of pp —?• jjj events with the same generator configuration as 
described previously. We also generated similar samples at y/s = 7 TeV for validation purpose : we 
compared the differential cross-sections we obtained to the ones measured by CMS . Our simple 
prediction is consistent with the reference (within 10%) at low invariant mass, but overshoots the 
reference by a factor up to 3 at large invariant mass. We didn’t correct our prediction to take 
into account this observation, which means that the results we present should be conservative as 
the level of background is clearly overestimated ®. Reducing the background by a factor 3 would 
only improve the limits on the effective coupling by about 30%, with essentially no impact on the 
qualitative conclusions set from these first projections. 

We also extended the invariant mass range probed, the upper bound increasing from 3270 to 
5540 GeV: we chose the latter so that the expected QCD background yield in the last bin for 
{^/s = 13 TeV, L = 5 fb“^) is 0(10) events, i.e. similar to the 7 TeV case. This extension is quite 
helpful as the sensitivity to the Z' signal comes mostly from the high end of the invariant mass 
spectrum. We used in the new range a variable-width binning (^^ = 6%), identical to the one 

® This only applies to the 13 TeV projections, which rely on the LO prediction, but not the reinterpretation of the 
7 TeV results which use the more accurate prediction from |23| . 
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used for the CMS measurement (Fig. 3(a)), which was chosen as to minimize the impact of the 
finite jet energy resolution on the measurement. 


We used a very similar setup to the one aforementioned to perform the hypothesis tests gauging 
the significance of a potential signal, only switching to a two-sided test statistic instead. The 
uncertainty on the background prediction is strongly inspired by the CMS results (cf Fig. 1 and 
3 in [23]), in which the theoretical uncertainties are comprised between 10 and 20%, to which 
should be added the experimental uncertainties dominated by sources related to the jet energy 
scale (JES, 5 to 30%). We decided to assign a flat uncertainty of 20%, counting on a future 
reduction of JES uncertainty for high pT jets solely based on increased statistics for performance 
measurements. The correlation of uncertainties between the different invariant mass bins is not 
straightforward though: if one assumes a fully correlated uncertainty, it might lead to an overly 
optimistic significance as the associate nuisance parameter can be strongly constrained in the MLE 
fit due to the bins at low invariant mass that have large statistics and are signal-free - an instance of 
so-called “profiling” which is undesired here. On the other hand, fully uncorrelated uncertainties 
may also lead to a too optimistic significance, as the signal generally spans several bins. We 
therefore adopted a conservative compromise, assuming full correlation but reducing the size of the 
uncertainty in the bins with low signal yields. Specifically, the uncertainty in each bin was set to 
20% X {Si/Bi)/ max(S'j/Sj), where Si and Bi are the respective signal and background yield in a 
particular bin i. Consequently, different uncertainty profiles are used for different Z' masses. 

Fig. I presents the discovery potential for integrated luminosities of 5 and 100 fb for two 
common levels of significance {3a or 5a) expressed in terms of Gaussian standard deviations. One 
can see that even with a luminosity not exceeding the one used for the CMS measurement at 7 
TeV, it would already be possible to probe free couplings up to one order of magnitude smaller 
than those excluded by the current searches. 

To finish, we’d like to mention existing searches |22l |29] for micro black holes or string balls, 
which select events with several jets and look at the Ht spectrum, the scalar sum of the jets px- 
While this signature is closely related to our scenario and could potentially be quite sensitive to 
Z' production, a quick estimate obtained from the measured Ht spectrum in [22] (Fig. 1) showed 
that these searches are less competitive than the results based on the three-jet invariant mass that 
we present here. 
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SM (NLO) 

Z' 300 GeV 

Z' 500 GeV 

Z' 800 GeV 

Z' 1.6 TeV 

Z' 3 TeV 

ys = 8 TeV 

~ 1.3 fb 

2.8 pb 

0.36 pb 

55 fb 

5.9 fb 

0.28 fb 

ys = 13 TeV 

9.2 fb [ni 

0.57 ^b 

74 pb 

11 pb 

1.2 pb 

57 fb 


TABLE III; Production cross-sections of the four-tops process pp —>■ tttt in Standard Model (leftmost column) 
and via Z' mediation (not including the SM contribution) with an effective coupling dg/M"^ = 10“®. 

V. SENSITIVITY IN THE FOUR-TOPS CHANNEL 


Measurements of rare Standard Model processes can be powerful tools to search for new physics 
in an indirect way. In our case, the Z' boson, through its coupling to gluons, might play an indirect 
role in QCD physics and particularly in top quark physics. We estimated the leading order Z'- 
mediated contributions to the cross-sections of a few SM processes involving top or bottom quarks, 
such as pp —)• ttg, pp —)• ttbb, pp —)> bbbb or pp —)• tttt. It turns out that potential Z' contributions 
manifest themselves most visibly in the latter process, the associate production of four top quarks, 
thanks to the very small corresponding SM cross-section (about 1 fb at ^/s = 8 TeV). This is 


illustrated in Table III, which provides cross-sections for a few Z' masses. The cross-sections were 
computed at leading order with MadGraph as described in the previous section. As it is especially 
relevant here, one should note that the cross-sections of the Z'-mediated contributions were seen not 
to depend on the Z' decay width, the resonant contributions being suppressed by the Landau-Yang 
theorem. Furthermore, we checked that interferences between Standard Model and Z'-mediated 
contributions are negligible (below 5%). 


The peculiar signature of four top quarks appears in various BSM scenarios such as Super- 
symetry or new heavy quark generations, and is looked for at the LHC |30II33j . To evaluate the 
current sensitivity of these searches to Z' bosons, we reinterpreted the results obtained in [31] in 
terms of upper limits on the Z' effective coupling to gluons, as was done in the previous section. 
In this publication by the CMS collaboration, events with four top quarks were looked for in the 
collision data produced by the LHC in 2012 {^/s = 8 TeV, L = 20 fb“^), using final states with one 
isolated lepton and jets. Minimal requirements were placed on the missing transverse momentum in 
the event (Li™®® > 30 GeV), the number of jets (at least 7 with pp > 30 GeV), and the scalar sum 
of selected jets and leptons transverse momentum {Hp > 400 GeV), after what a boosted decision 
tree (BDT) was used to discriminate signal from background events, in four distinct categories 
(electron or muon, exactly 7 or > 8 jets). The hypothetical signal yield was extracted through a 
combined fit of the final BDT discriminant distributions in the four categories. 
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FIG. 5: Discovery potential at LHC during the Run-2 phase, in terms of the Z'-gluon coupling, together 
with current constraints established here and in m- 

The absence of excess over the expected background lead to the establishment of an upper 
limit on the cross-section of the pp —)• tttt process of 32 fb, that is 24 times the signal strength 
of the Standard Model process. Assuming that kinematic distributions do not vary significantly 
between Standard Model and Z'-mediated contributions, this limit can be directly translated into 
a limit on the Z'-gluon coupling. This new constraint is represented in Fig together with the 
limits obtained from the three-jets final state and the dark matter-related constraints. One can 
notice that the four tops final state brings a useful complementarity to the other channels for very 
heavy Z' (above 3 TeV), since the upper limit on the effective coupling increases only linearly 
with the Z' mass, while the sensitivity in the three-jets channel vanishes quickly when the Z' mass 
approaches the collider center-of-mass energy. 


Projections for the next LHC run at ^/s = 13TeV 

For this channel as well, we estimated the discovery potential that may be achieved during the 
LHC second phase of exploitation at ^/s = 13 TeV. Heavy final states, such as those with four top 
quarks, will largely benefit from the increased center-of-mass energy. As can be seen in Table 0 
this is true not only for the Standard Model process (cross-section increased by a factor 7), but 
also for Z'-mediated contributions in much greater proportions (a factor 200 increase). 

Taking the simplest but most approximate approach to extrapolate the current experimental 
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results, we relied on the distributions of the BDT discriminant for the main backgrounds 

(tt+light jets, tt + cc/bb) and Standard Model tttt contribution that are presented in the CMS 
publication (Fig. 3 in [31]). We reweighted these distributions by the ratios of the leading order 
cross-sections between ^/s = 13 and 8 TeV, which are respectively ~ 5.3 {tt+ light jets), ~ 7 
{tt -|- ccjbb) and ~ 7.4 {tttt). We then used these reweighted distributions to evaluate the signal 
significance at ^/s = 13 TeV through a combined fit of the signal strength in the four categories, the 
hypothesis test being performed again with the HistFitter software. We assigned global systematic 
uncertainties of 20% on the tt+jets, and 50% on the tt + cc/bb background predictions, reflecting 
the total uncertainties mentioned in m- We first checked that this configuration allowed us to 
reproduce the CMS analysis result at 8 TeV : we indeed obtained an upper limit on the pp —)■ tttt 
cross-section only differing by 5% from the reference. 

Fig. in presents the estimated sensitivity at 13 TeV, again for integrated luminosities of 5 and 
100 fb“^. The sensitivity in terms of the effective coupling is improved by a factor 2-3 with respect 
to the sensitivity in this channel at 8 TeV, but does not allow to go beyond the upper limits already 
set by the 7 TeV experimental results in the three-jets channel apart at large Z' masses (above 2 
TeV), a region the four-tops channel is the best tool to probe. 


VI. CONCLUSIONS 

In this paper, we studied the effect of the presence of a gluophilic Z' on multi-tops event at the 
LHC, and in particular three jets and four tops events. The only coupling parameter of the model 
- the effective coupling involved in dimension six operators - has been constrained in detailed from 
experimental searches released at 7 and 8 TeV by CMS along these two different channels. The 
latter restrictions led us to evaluate for which values of the parameter space a gluophilic Z' could 
be discovered during the next run of the LHC. In particular, the analysis of the three jets invariant 
mass could provide a clear signal (> 5a) for masses of the Z' above 300 GeV. In addition, the four 
tops events analysis could furnish a potential of discovery for heavy Z' (above ~ 2 TeV). Existing 
constraints coming from the study of dark matter |10j where the studied Z' represents a possible 
mediator of the latter were furthermore added to the analysis, showing a tension with potential of 
discovery for low masses of the Z' : masses lower than 300 GeV are disfavored to be detected from 
this perspective. 

In conclusion, it is important to notice that a possible combination of both measured excesses 
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in four top channels and three jets analysis (for large invariant masses) during the next run of the 
LHC could provide a clear signal of the presence of a gluophilic Z'. 
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